Resistively-shunted superconducting quantum point contacts 
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We have studied the Josephson dynamics of resistively-shunted ballistic superconducting quantum 
point contacts at finite temperatures and arbitrary number of conducting modes. Compared to 
the classical Josephson dynamics of tunnel junctions, dynamics of quantum point contacts exhibits 
several new features associated with temporal fluctuations of the Josephson potential caused by 
fluctuations in the occupation of the current-carrying Andreev levels. 
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Andreev levels with energies below the superconduct- 
ing energy gap A are known to carry the stationary 
supercurrent in ballistic superconducting point contacts 
y-g|, and are also responsible for the dc and ac current 
flow at small bias voltages V <C A/e [§|-§|. Occupa- 
tion factors of these states vary randomly, either due to 
the thermal fluctuations in the contact electrodes [0,||, 
or due to the probabilistic nature of the Landau-Zencr 
transitions [^ induced by the bias voltage in contacts 
with finite reflection coefficients |Q]. These random fluc- 
tuations in the occupation factors of the sub-gap states 
lead to a large supercurrent noise. At finite bias volt- 
ages, the noise can be interpreted as the shot noise of the 
large charge quanta of magnitude 2A/T^ which is equal 
to the charge transferred through the contact during one 
period of the Josephson oscillations. Recently, this noise 
has been observed in experiments with high-transparency 
tunnel junctions pO[ . 

In a point contact with an external environment of 
finite impedance, the supercurrent noise leads to fluc- 
tuations of the voltage across the contact, and affects 
the dynamics of the Josephson phase difference tp. In 
the case of quantum point contacts with few propagating 
electron modes, for example those fabricated with the 
controllable break junction technique |ll|-[l^], the typ- 
ical external impedance is much smaller than the con- 
tact resistance. In this case the effects of the external 
impedance are signiflcant only at small bias voltages, 
where dynamics of ip can be described within the adi- 
abatic approach Q. The aim of this work is to study 
adiabatic flnite-temperature dynamics of superconduct- 
ing quantum point contacts, taking into account the su- 
percurrent noise caused by fluctuations in the occupation 
factors of the Andreev levels. 

We consider the standard model of a point contact, 
characterized as a short ballistic constriction support- 
ing N electron modes propagating between two identical 
bulk superconductors. At small bias voltages, V <C A/e, 
the current through such a constriction is carried by 
the two levels (per each electron mode) with energies 
±e((^), where e{ip) — Acos((/j/2). We assume that 
the external impedance of the contact reduces to the 
frequency-independent constant R in the relevant fre- 



quency range around the typical Josephson oscillation 
frequency of 2eIcR/fi, where Ic = NAe/h denotes the 
zero-temperature critical current of the contact. 

The basic result of the adiabatic theory [^ of dynamics 
of superconducting point contacts is that each occupied 
Andreev state contributes energy e(</j) or —e{ip) to the 
Josephson coupling energy of the point contact, while the 
occupation probability w of the state is determined by 
the quasiparticlc transitions between this state and the 
electrodes of the contact. These transitions are described 
by the simple rate equation: 



w = 7(e) [/(e) -w] 



(1) 



where /(e) is the equilibrium Fermi distribution, e = 
e{ip(t)), and the transition rate 7(e) depends on the in- 
elastic or pairbreaking effects which give rise to the finite 
subgap density of states in the electrodes H . Equation 
Imj shows that the role of quasiparticlc transitions is to 
drive the occupation probabilities of Andreev states to- 
wards the equilibrium. 

The time evolution of the Josephson phase differ- 
ence (p{t) affects the occupation probabilities of Andreev 
states through the phase dependence of their energies 
±e{ip). At finite temperatures the phase dynamics is dif- 
fusion along the Josephson potential driven by the ther- 
mal noise generated by the external resistance R. Such 
a diffusion is governed by the following set of Fokker- 
Planck equations for the probability density a: 
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where k = —N,—N+l,...,Nis the difference of the 
numbers of the occupied Andreev levels that carry pos- 
itive and negative currents. It characterizes different re- 
alizations of the Josephson potential, U^''' = —ke{p) — 
(hlQ/2e)p, where /q is the external bias current (see Fig. 
1). The relaxation term r{(Tfe} in eq. (g) describes the 
transitions from one realization of Josephson potential to 
another, induced by the quasiparticlc transitions (|l|) be- 
tween the Andreev states and the electrodes of the point 
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FIG. 1. Equivalent circuit of the resistively shunted su- 
perconducting quantum point contact, and three reahzations 
of the Josephson potential for a single-mode ballistic contact. 
The phase diffusion along this potential is characterized by the 
three-component probability density ak{<p) which should be 
matched continuously between the points ifi — and (p = 27r 
- see eqs. (|9|) and (|ic|). 

contact. Taking a sum over transitions from, and into, 
all of 2N Andreev levels we get for F: 

T{ak} = 7(e)/(e) [{N + k + l)afe+i - {N + k)ak] + 

7(e)/(-e) [{N-k + l)afc_i - {N - k)ak] . (3) 

Since the dependence of all "internal" properties of 
the point contact (for instance, Andreev state energies 
±e((/j)) on (p is periodic with the period 27r, it is con- 
venient to restrict (p to the interval [0, 27r]. The diffu- 
sion equation (0) should then be complemented with the 
boundary conditions for ct's at the ends of this interval, 
if = 0, 27r. The boundary conditions are somewhat com- 
plicated by the fact that when the phase ip approaches 
or 2n the energies of the Andreev state hit the gap 
edges ±A. As a result, the quasiparticle transitions be- 
tween these states and electrodes become "absolutely" 
efficient, leading to instantaneous equilibration of the oc- 
cupation probabilities of Andreev states . (More accu- 
rately, one should say that this equilibration process has a 
time scale that is negligible in comparison to the period 
of the Josephson oscillations.) Since different Andreev 
states are occupied independently one from another, this 
implies for the probability densities: 



a,ip = 0) = A[fi-A)f+'^[fiA)]^-'^C^+'' , 



]N~k^N+k 



au{^ = 2n) = A'[!{A)r+'[f{-A)f~'a 



2N 



(4) 
(5) 



where A, A' are some constants, and C,™ = nl/[m\{n — 
m)!]. 

Relations (Q) and ^ can serve as the boimdary con- 
ditions for the diffusion equations (0). One more condi- 
tion follows from the requirement that the components 
akif) of the probability density diffusing on the different 
branches U^''\ip) of the Josephson potential be contin- 
uous. Taking into account that the states of the point 
contact with (p and ip + 2tt should be identical, we can 
write this condition as A' = A. Combining this result 
with eqs. (^) and (||), we see that one of this equations 
can be replaced with the following boundary condition: 



ak{(p = 27r) = CT_fc((y9 = 0) 



(6) 



Equations (Q) and (^), together with the usual normal- 
ization condition 



V / dipak{(p) = 1, 

I, "'0 



(7) 



determine completely the solution of the diffusion equa- 
tion (g), and thus allow us to find all the characteristics 
of the point contact. Particularly, in the steady-state 
case, eq. (||) reduces to the equation djk/dip = Fjcrfc}, 
and can be solved explicitly to obtain the dc voltage V 
across the point contact: 



v = — Vjfc. 

e -^-^ 

k 



(8) 



Equations (§) ~ (o) show that the time evolution of the 
ip depends, in a non-trivial way, on the number N of the 
propagating electron modes in the point contact. Below 
we analyze two opposite limits of such a dependence, N = 
1 and N '3> 1- We also assume that the subgap density 
of states in the electrodes of the point contact is small 
and we can neglect the quasiparticle exchange term F in 
eq. (0). We start with the case of a single-mode quantum 
point contact with A^ = 1 . We have then three diffusion 
equations (H) with fc = 0; ±1 (see Fig. 1). The boundary 
conditions (||) and (^) are: 



ai(0) 
and 



-ao(0)e^/^, a_i(0) = -ao(0)e 



-A/T 



f7±i(27r) = f7^i(0), ao(27r) = (To(0). 



(9) 



(10) 



The stationary version of eq. (0), jk{f) — const, can 
be integrated directly, as in the case of classical Joseph- 
son junctions [|l4|Jl5|]. T aking into account the boundary 
conditions (H) and (|lO|) we find then the cr's: 



CT±i('/5) = ^exp{- 
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T 
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where K± = J^^ dipexp{U^'^^\ip)/T}, and (Tq is a con- 
stant density ctq (ip) that is determined by the normahza- 
tion condition 
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(12) 



In terms of ctq , the dc voltage V across the point contact 
is: 



V = 2TraoIoR 



i + H'^ + T^'" --"""'"' 



(13) 



Figure 2 shows the current I flowing through the point 
contact, I = Io — V/R (relation between I and Iq is illus- 
trated by the equivalent circuit in Fig. 1), as a function of 
the voltage V calculated from eqs. (|l])-(ll^) for several 
temperatures. The zero-temperature I — V characteristic 
(the upper curve in Fig. 2) can be represented analyti- 
cally by the following dependence of the voltage V on the 
bias current Iq M: 



V 



vri?(/2-/2)i/2 



4 arctan y/{Io + Ic)/{Iq ~ h) 



h > Ic 



(14) 



where Ic = elS.j'h is the zero-temperature supercurrent of 
the point contact. This equation describes the transition 
from the supercurrent at 1^ = to a constant current 
2Ic/tt at V^ :^ RIc. The large- voltage dc current arises 
from the strongly non-equilibrium occupation of the An- 
dreev levels, and is a hallmark of multiple Andreev re- 
flections in a ballistic point contact. 

Both the zero-voltage supercurrent and the large- 
voltage current are suppressed by temperature. For in- 
stance, in the limit of large temperatures, T 3> A, we get 
from eqs. (0)-(|ll): 
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(15) 



This equation describes well the lowest curve in Fig. 2. 
Transition from low to high temperatures can be traced 
analytically at large bias currents, Iq ^ /c, when eqs. 
(|l|)-(|| give: 



2 A 

/= _/^tanh(— ), 

IT Zl 



(16) 



in agreement with Fig. 1 and large-current limit of eq. 

At small bias currents Iq < Ic, and small temperatures 
T <IC A, the point contact is in the supercurrent state 
and the bias current flows almost completely through 
the contact, I ~ Iq. The voltage across the contact is 
exponentially suppressed and is associated with the rare 




FIG. 2. DC I — V characteristics of the single-mode super- 
conducting ballistic quantum point contact at several temper- 
atures. All components of the current are suppressed when 
temperature T increases above the energy gap A. The vari- 
ous limits of these curves are given analytically by eqs. (114) 
-(111). 



"phase-slip" events (^ — > (/? -t- 27r. In this regime, and for 
bias currents on the order of the critical current Ic, we 
get from eqs. (0)-(0|): 



^=^(l-*Y/^cxp{-^[(l-z2)V^ 



—i{ — 



^W)]}, 



' = z- 



(17) 



This equation corresponds to the classical thermal ac- 
tivation over the maximum of the Josephson potential, 
similar to the one found in classical Josephson junctions 

Now we consider point contacts with large number of 
transverse modes, TV ^ 1. The probability density ak 
is strongly peaked in this case as a function of k at its 
average value: k — (2i(; — 1)N, where w is the occupa- 
tion probability of the Andreev states carrying positive 
current. This means that we can neglect fluctuations of 
the Josephson potential around the average potential: 



U{ip) ^Ne{ip)il-2w) 



(18) 



The probability w, which determines the shape of the po- 
tential (iq), depends sensitively on the equilibration rate 
7 in the rate equation (H^). We start by considering the 
case of an ideal BCS superconductor in which 7 is com- 
pletely suppressed in the relevant sub-gap energy range, 
and equilibration of the Andreev states occurs only at 



(f — 0, 2tt. In this case the probabiUty w is equal to 
/(—A) or /(A) depending on whether the phase diffuses 
into the interval [0, 2tt] through one {(p = 0), or the other 
(if = 27r), end of this interval. This means that although 
we neglect fluctuations of the Josephson potential, we 
still have two branches of this potential due to the non- 
equilibrium occupation of Andreev states. The potential 
is Ej cos{(p/2) on one branch, and —Ej cos{ip/2) on the 
other, where Ej = NAtauh{A/2T). The phase diffu- 
sion along these two branches can be described with the 
two equations in (0) with k = ±1, if we replace the en- 
ergy gap A in U^^'{ip) with Ej. However, the boundary 
conditions for a±i are now quite different from those we 
used for the single-mode junction - see eqs. (§) and (10). 
As we discussed above, the phase diffuses along one or 
the other branch of the potential depending on whether 
it enters the interval [0, 2tt] through one or the other end. 
This dynamics is accounted for by the following bound- 
ary conditions: 



ai(27r)=a_i(0) = 0, a_i(2^) = ai(0) . (19) 

Equations (0) with the boundary conditions (O) give: 
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a_i(¥.) = ^e-"°/-^exp{- - ^ ^"^ } 

X / Vexp{ -^}, (21) 







T 



and 



V = ai 



(22) 



where ct is a constant that is determined by the normal- 
ization condition (M). 

Figure 3 shows the I — V characteristics of the 
multi-mode point contact calculated from eqs. (pQ)- 
(E3). The zero-temperature curve is given by the same 
eq. ([^ as for iV = 1, but with /^(T) = eEjjh = 
iV(eA/S,)tanh(A/2r). (Note that in the multi-mode 
contact, Ej is much larger than A, and the "zero- 
tempertature" I — V characteristics with T ^ Ej can 
still correspond to T ::^ A.) Since the current in Fig. 3 is 
normalized to the tempcraturc-dcpcndcnt critical current 
/c, the large- voltage limit of the current in Fig. 3 is in- 
dependent of temperature and equal to 2/c/7r in contrast 
to the single-mode case shown in Fig. 2. 

At large temperatures, T 3> £',/, the current / can be 
found analytically from eqs. (EQ)-(E2h: 



/ = 



2E 



TT 1 



^coth( — ), 1 = ^7^. 



eT 



(23) 



This equation agrees with the lowest curve in Fig. 3. At 
low temperatures, T ^ Ej^ and small bias currents /q < 




FIG. 3. DC I — V characteristics of the multi-mode quan- 
tum point contact between two BCS superconductors at sev- 
eral temperatures. For discussion see text. 



/c, the contact is in the supercurrent state and the voltage 
is generated only by the rare thermal activation over the 
Josephson energy barrier as in the single mode case. For 
not-too-small bias currents, the voltage is given by the 
same eq. dl^) with A replaced by Ej. At very small 
currents, IqJIc < {T/EjY/'^, the voltage is: 



T:hln 



y = 2/,i?e-2S'/^sinh(^) 



eT 



(24) 



Finally, we qualitatively discuss the effect of the fi- 
nite rate 7 of quasiparticle exchange between the An- 
dreev states and contact electrodes. At small voltages 
V <C Ti'^ I e the phase evolves slowly and the quasiparticle 
transitions maintain equilibrium occupation of Andreev 
states. This means that the Josephson potential (|l^) re- 
duces to 



U{lp) = -Ne{ip) tanh( 



2T ' 



(25) 



The main result of such a reduction is that now the 
potential depends only on the instantaneous value of 1^, 
and therefore has only one branch instead of two as in 
absence of the quasiparticle relaxation. The exact shape 
of the potential (Eq) depends on temperature, changing 
from -iVA | cos((p/2) | at T < A to -(iVAV4T) cos(p 
at T ^ A. In this temperature range, T ~ A, the tem- 
perature is still negligible on the scale of the Josephson 
potential, and we can use the zero-temperature version 
of eqs. (0) to describe the dynamics of Lp. Combining 
this equation with the boundary condition (T(27r) = cr(0) 



appropriate for the situation with one branch of the po- 
tential, we get for the voltage V: 



V ^2ttR 



27r 



dip 



Io-Is{f) 



(26) 



where Ig = {2e/h)dU{ip)/dip. In the two hmits T > A 
and T ^ A, eq. (|2q ) gives the standard result of the 
resistively-shunted-junction (RSJ) model, V — R{Iq — 
^c)^ J where Ic =niaxy[/s(^]. It can be shown by nu- 
merical integration of eq. (Eq) that deviations from this 
result are smaller than a few percent for arbitrary A/T 
ratio. 

When the temperature becomes nonvanishing on the 
scale of Josephson potential, it is already much larger 
than A and the potential (ESh becomes the regular 
Josephson potential —Ej cosip. This means that the 
quasi-equilibrium dynamics of the multi-mode ballistic 
contacts can be described with the classical RSJ model 
M,15[. When the voltage V across the contact is suffi- 
ciently large, V 3> h^/e, this quasi-equilibrium dynamics 
goes over into the non-equilibrium one, described by Fig. 
3 and eqs. (0), (||). 

In summary, we have presented a new approach to 
the description of the classical Josephson dynamics of 
superconducting quantum point contacts with arbitrary 
number of conducting quantum modes in a resistive en- 
vironment. New characteristic features of the Joseph- 
son dynamics in these contacts include fluctuations of 
the Josephson potential caused by fluctuations in the 
occupation of the current-carrying Andreev levels, and 
low-voltage excess current associated with the non- 
equilibrium occupation of these levels. 
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